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ABSTRACT 

We present new VLA C+D-array H I observations and optical and NIR imaging of the well known 
interacting system NGC 4038/9, "The Antennae". At low spatial resolution (~ 40"), the radio data 
reach a limiting column density of ^ 10^^ cm^^ (2.5(t), providing significantly deeper mapping of the 
tidal features than afforded by earlier Hi maps. At relatively high spatial resolution 10"), the radio 
data reveal a wealth of gaseous sub-structure both within the main bodies of the galaxies and along 
the tidal tails. In agreement with previous H I studies, we find that the northern tail has H I along its 
outer length, but none along its base. We suggest that the H I at the base of this tail has been ionized 
by massive stars in the disk of NGC 4038. The gas in the southern tail has a bifurcated structure, 
with one filament lying along the optical tail and another running parallel to it but with no optical 
counterpart. The two filaments join just before the location of several star forming regions near the end 
of the tail. The H l velocity field at the end of the tail is dominated by strong velocity gradients which 
suggest that at this location the tail is bending away from us. We delineate and examine two regions 
within the tail previously identified as possible sites of a so-called "tidal dwarf galaxy" condensing out 
of the expanding tidal material. The tail velocity gradients mask any clear kinematic signature of a 
self-gravitating condensation in this region. A dynamical analysis suggest that there is not enough mass 
in gas alone for either of these regions to be self-gravitating. Conversely, if they are bound they require 
a significant contribution to their dynamical mass from evolved stars or dark matter. Even if there are 
no distinct dynamical tidal entities, it is clear that there is a unique concentration of gas, stars and 
star forming regions within the southern tail: the H I channel maps show clear evidence for a significant 
condensation near the tail star forming regions, with the single-channel H I column densities higher than 
anywhere else in the system, including within the main disks. Finally, the data reveal H I emission 
associated with the edge-on "superthin" Scd galaxy ESO 572-G045 that lies just beyond the southern 
tidal tail of The Antennae, showing it to be a companion system. 

Subject headings: galaxies: evolution — galaxies: individual (NGC 4038/9) — galaxies: interactions 
— galaxies: ISM — galaxies: kinematics and dynamics — galaxies: peculiar 



1. INTRODUCTION 

NGC 4038/9 (= Arp 244), aptly nicknamed "The An- 
tennae" , is one of the best-studied examples a galactic col- 
lision. The long tails which distinguish this system are em- 
blematic of violent tidal interactions between disk galax- 
ies of similar mass (Toomre & Toomre 1972, Schweizer 
1978). Tidal tails originate from the outer regions of galac- 
tic disks and are often rich in neutral atomic hydrogen. 
These features develop kinematically and their velocity 
fields bear the imprint of the encounter process, so they 



provide crucial information for dynamical modeling of in- 
teracting galaxies. We therefore targeted NGC 4038/9 for 
high-spatial and velocity resolution H I observations (full- 
width at half maximum, FWHM, as fine as ~ 10"; Aw ^ 
5 km s~^) in order to obtain the data needed to constrain 
future numerical simulations of this system. 

NGC 4038/9 was the first major merger to be mapped 
in the 21-cm line of neutral hydrogen. The presence of 
significant quantities of neutral hydrogen within both the 
main bodies and the tails of NGC 4038/9 was established 
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from single-dish measurements (Peterson & Shostak 1974, 
Huchtmeier & Bohnenstengel 1975). Subsequently, NGC 
4038/9 has been the subject of Hi synthesis mapping 
observations using the Westerbork Synthesis Radio Tele- 
scope (WSRT; van dcr Hulst 1979a), the Very Large Ar- 
ray (VLA; Mahoney, Burke & van der Hulst 1987) and 
the Australian Telescope Compact Array (ATCA; Gordon, 
Koribalski & Jones 2001). The observations reported in 
this paper provide a more detailed mapping of the tidal 
features than afforded by previously published H i maps, 
representing an improvement in spatial and velocity reso- 
lution by over a factor of two. 

Other spectral line studies have mapped the kinematics 
of NGC 4038/9 using CO spectral hues (Stanford et al. 
1990, Wilson et al. 2000, Gao et al. 2001, Zhu 2001) and 
optical emission lines (Burbidge & Burbidge 1966, Ruben, 
Ford & D'Odorico 1970, Amram et al. 1992). However, the 
molecular and ionized gas is confined to the inner disk re- 
gions, and their kinematics do not constrain the encounter 
geometry as well as that of the tidal tails. 

At a heliocentric velocity^ of 1630 km s~^ (correspond- 
ing to a velocity relative to the local group of 1440 km 
s""'^). The Antennae is one of the nearest on-going major 
mergers. To be consistent with the majority of the recent 
work on NGC 4038/9, we derive physical quantities by 
adopting a Hubble Constant of 75 km Mpc~^, yield- 
ing a distance of 19.2 Mpc. However, a recent analysis of 
HST VI photometry of individual stars within the south- 
ern tidal tail of NGC 4038/9 by Saviane, Rich & Hibbard 
(2001) suggests a much smaller distance of 13.8±0.5 Mpc. 
This distance is derived by identifying luminous red stars 
within the tail with the red giant branch of an old, metal 
poor population. There are some caveats to this deriva- 
tion, for which we refer the reader to Saviane et al., so 
for the purposes of the present work, we will use the more 
commonly accepted value of 19.2 Mpc. 

At this distance, the tidal tails of NGC 4038/9 extend 
some 65 kpc in radius and measure ~ 110 kpc from end to 
end. The bodies of the galaxies are sites of extensive star 
formation, producing an IR luminosity of log(LiR,) = 10.76 
and an IR to blue luminosity ratio of Ljr/Lb — 6. Long- 
wavelength studies from the NIR to Radio suggest that 
this luminosity is powered by an active system-wide star- 
burst, with no indication of a significant contribution by 
an AGN (Hummel & van der Hulst 1986, Vigroux et al. 
1996, Fischer et al. 1996, Kunze et al. 1996, Mirabel et al. 
1998, Nikola, T. et al. 1998, Laurent et al. 2000, Haas 
et al. 2000, Xu et al. 2000, Neff & Ulvestad 2000, Mengel 
et al. 2001). While the IR luminosity is an order of magni- 
tude lower than the most active star forming mergers (the 
so-called Ultraluminous Infrared Galaxies), it is still a fac- 
tor of ~5 higher than non-interacting galaxy pairs (e.g. 
Kennicutt et al. 1987, Bushouse 1987, Bushouse, Lamb & 
Werner 1988), with an inferred star formation rate (SFR) 
of ~ 20AfQyr-i (Evans, Harper & Helou 1997; Zhang, 
Fall & Whitmore 2001). 

NGC 4038/9 is the nearest system with an identified 
dwarf-galaxy sized concentration of gas, light and young 
stars embedded within a luminous tidal tail (Schweizer 
1978, hereafter S78; Mirabel, Dottori, & Lutz 1992, here- 
after MDL92). Such concentrations have been known for 

^Heliocentric velocities are quoted throughout this paper. 



quite some time (e.g. Zwicky 1956), but have only re- 
cently received detailed attention. Several numerical stud- 
ies lend support to the hypothesis that such systems are 
self-gravitating and may evolve into independent dwarf 
galaxies (Barnes & Hernquist 1992, Elmegreen, Kaufmann 
& Thomasson 1993), but the supporting observational 
evidence is primarily circumstantial. These observations 
show that the optical condensations within tidal tails con- 
tain many young stars (S78, MDL92, Hunsberger, Charl- 
ton, & Zaritsky 1998, Weilbacher et al. 2000, Iglesias- 
Paramo & Vi'lchez 2001, Saviane, Rich & Hibbard 2001) 
and have global properties, such as size, luminosity, H I 
mass. Hi velocity dispersion, and CO content, in common 
with dwarf Irregular galaxies (Mirabel et al. 1992; Hibbard 
et al. 1994; Smith & Higdon 1994; Due & Mirabel 1994, 
Due et al. 1997, 2000; Smith & Struck 2001; Braine et al. 
2000, 2001). Recently, some tidal dwarfs have even been 
found to harbor "Super Star Clusters" (Knierman et al. 
2001). For these reasons, such concentrations are often 
referred to as "Tidal Dwarf Galaxies" (hereafter TDGs). 
However, strong kinematic evidence that these concentra- 
tions are indeed self-gravitating is lacking, and their evolu- 
tion into dwarfs is therefore still in question (see Hibbard, 
Barnes, & van der Hulst 2001, hereafter Paper II, for more 
details). Another objective of the present observations is 
therefore to search for the kinematic signature of a self- 
gravitating mass in the vicinity of the putative TDG. 

The paper is organized in the following manner. In §2 
we present details on the observations and data reduction 
procedures. In §3 we describe in turn the morphology and 
kinematics of the tidal tails, sub-structure within the tails, 
the region in the vicinity of the TDG, the inner disks, and 
the companion galaxy. In §4 we try to explain various as- 
pects of the H I tidal morphology, and look for evidence for 
a kinematically distinct entity at the end of the southern 
tail. Our conclusions are presented in §5. 

2. OBSERVATIONS 
2.1. VLA Hi Observations 

The Antennae was observed in two separate array con- 
figuration of the VLA radio interferometer: the most com- 
pact configuration (1 km D-array) in June of 1996, and the 
hybrid CnB-array in June of 1997. For the hybrid array, 
the antennae along the E and W arms of the array are at 
the standard C-array (3 km) stations while those of the 
N arm are at the more extended B-array (9 km) stations. 
This compensates for the foreshortening of baselines to the 
north arm antennae that occurs for sources at low decli- 
nation. The short baselines of the D-array provide the 
highest sensitivity to extended gas, while the longer base- 
lines of the CnB-array allow resolutions as fine as ~ 10". 
NGC 4038/9 was observed for a total of 3.5 hours in the 
D-array and 8 hours in the CnB-array The details of these 
observations are tabulated in Table 0. 

We observed using the 21-cm spectral line mode with 
the system tuned to a central frequency corresponding to 
a heliocentric velocity of 1630 km s^^, the central velocity 
of the existing H i synthesis observations (Mahoney et al. 
1987). Since the atomic gas within tidal tails typically 
has narrow linewidths (tT„ ~ 5 — 10 km s~^; Hibbard et al. 
1994, Hibbard & van Gorkom 1996), the correlator mode 
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was chosen to give the smallest channel spacing while still 
covering the velocity spread of existing H i observations 
(1360-1860 km s"^; van der Hulst 1979a, Mahoney et al. 
1987). New broader-bandwidth Hi (Gordon et al. 2001) 
and CO (Gao et al. 2001) observations reveal that there 
is both atomic and molecular gas over the velocity range 
from 1340-1945 km s^^. We discuss the implications of 
this result below. 

Using on-line Hanning smoothing, the chosen set-up 
provides 127 independent velocity channels over the 3.125 
MHz bandwidth, resulting in a channel spacing of 5.21 km 
s^^. The pointing center was chosen to place the ends of 
the optical tails at equal distances from the phase center. 
This position is 2.'5 south of the main bodies, and places 
the ends of the optical tails at a radius of ^ 10', which 
is approximately the 75% sensitivity point of the primary 
beam (see Figure 0). 

The data were calibrated, mapped, and "Cleaned" using 
standard techniques and procedures in the Astronomical 
Image Processing System (AIPS) (see e.g. Rupen 1999). 
The data from the separate arrays were calibrated and 
bandpass corrected independently, after which they were 
combined in the UV plane to produce the C-fD dataset. 
This latter dataset is used exclusively throughout this pa- 
per. The final calibration of this dataset was achieved by 
three iterations of a phase-only self-calibration. The re- 
sulting continuum map is very similar to others that have 
appeared in the literature (e.g. Hummel & van der Hulst 
1986; Neff & Ulvestad 2000, Gordon et al. 2001), and is 
not shown here. 

The continuum was subtracted from the velocity cube 
by fitting a first-order polynomial to the visibilities from 
the line free channels on either end of the bandpass. This 
is done in an iterative manner in order to determine which 
channels are free of line emission. The resulting fit was 
made using channels 11-20 and 104-111 of the 127 chan- 
nel cube, corresponding to the velocities 1859-1911 km s^^ 
and 1386-1423 km s^^. As mentioned above, cold gas has 
been found within both of these velocity ranges by Gor- 
don et al. (2001) and Gao et al. (2001). The gas at these 
velocities is of low surface brightness and totally confined 
to either the disk of NGC 4039 or the disk overlap region. 
Therefore, using this continuum range will not affect any 
of the measured properties of the tidal tails, which is our 
primary interest in this paper. However, the H i fiux we 
measure for the central region will underestimate the true 
flux, and the kinematics of this region will not be accu- 
rately mapped (although the fact that the line kinematics 
are weighted by intensity will somewhat mitigate the ef- 
fects of the missing velocity information). 

The continuum subtracted data were mapped using sev- 
eral values of the Robust weighting parameter, R (Briggs 
1995). This parameter can be varied to emphasize either 
the outer (i? < 0) or the inner {R > 0) regions of the 
UV plane (i.e., either small or large spatial scales). In 
practical terms, a lower R gives a finer spatial resolution 
at the expense of surface brightness sensitivity, while a 
larger R gives a higher surface brightness sensitivity at 
a slightly coarser spatial resolution. A high resolution 
data cube was made with R = —1, providing a surface 
brightness rms sensitivity of 1.3 mJy beam~^ at a reso- 
lution of 0FWHM = 11.4" X 7.4". This corresponds to a 
single channel column density limit (2.5ct per channel) of 



2.2xl0^"cm^^ averaged over the beam width of 1.1x0.7 
kpc^. A more sensitive intermediate resolution data cube, 
made with R = -1-1, provides a resolution of 21" x 15" 
(1.9x 1.4 kpc^) and a column density limit of 4x 10^^ cm^^. 
To further increase sensitivity to extended low column den- 
sity gas, a low resolution data cube was made by convolv- 
ing the R = +1 data cube to a resolution of 40" (3.7 
kpc), reaching a detection limit of 1.2x 10^^ cm~^. These 
data will be referred to in the following as the high, in- 
termediate, and low resolution data, respectively. These 
parameters are summarized in Table |^. The resulting per 
channel noise is close to theoretical once the higher sys- 
tem temperature due to the low declination of the source 
is taken into consideration. 

Visualization of the resulting data cubes and inter- 
comparison with the optical and NIR data was greatly 
facilitated by using the Karma visualization package 
(Gooch 1995). Besides allowing multiple images to be 
interactively compared, this package also allows a full 3- 
dimensional rendering of datacubes, which is a particularly 
powerful tool for disturbed systems. 

After a careful viewing, the data cubes were integrated 
over the velocity axis to produce the moment maps. In or- 
der to suppress the effects of noise, only data above a fixed 
threshold are used in the moment summation (the "cut- 
off technique" described by Bosma, 1981) using the AIPS 
task MOMNT. This task applies a user specified thresh- 
old to a version of the data cube smoothed in both space 
and velocity. If a pixel in the smoothed data cube passes 
the threshold, the corresponding pixel in the un-smoothed 
data is included in the moment analysis. This proce- 
dure favors low level emission that is extended in velocity 
and/or space over emission which is isolated. The output 
from this analysis is an integrated intensity map (zeroth 
moment), an intensity- weighted velocity map (first mo- 
ment), and an intensity- weighted velocity dispersion map 
(second moment). It should be noted that the first and 
second moment maps give an accurate representation of 
the mean H i velocity and line-of-sight velocity dispersion 
at a location only if the line profiles are single-peaked. The 
final zeroth moment maps were corrected for the primary 
beam attenuation before measuring Hi fluxes. 

We measure H i fluxes from an integrated intensity map 
constructed to match the new ATCA H i observations of 
Gordon et al. (2001), since those authors were able to do 
a proper line-free continuum subtraction. For the moment 
analysis, this involved applying a threshold of 2.1 mJy 
beam""^ to the data after smoothing with a 3x3 pixel box- 
car spatial filter and a five-channel Hanning filter. Our re- 
sulting fiux measurements are given in Table ^ The errors 
attached to the flux measurements have been determined 
from the single channel noise by properly taking into ac- 
count the number of independent channels and beams over 
which the flux was integrated. It should be noted that the 
error calculation does not include systematic errors, such 
as those introduced by poor continuum subtraction or in- 
complete cleaning. 
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2.2. Optical and NIR Observations 

The supporting optical observations were obtained dur- 
ing two nights of broad-band BVR imaging on the CTIO 
0.9m with the Tek512 CCD in April of 1991 and one night 
of near-Infrared (NIR) iC'-band imaging on the UH 88" 
with the QUIRC detector in January of 1995. The details 
of these runs are given in Table |^. For the optical run the 
conditions were photometric and the seeing was ^ 1". The 
f/13.5 re-imaging optics were used, giving a plate scale of 
0.44" pixel^^ and a field of view of 3.8'. The data were 
calibrated via observations of standards in the Graham E- 
regions (Graham 1982) observed on the same nights, with 
zeropoint errors {la) of 0.04™ in B and 0.02™ in both V 
and R. For the NIR run, the data were obtained through 
a K' filter (A = 2.11^m, AA = 0.35^m; Wainscoat & 
Cowie 1992; hereafter referred to simply as K). The f/10 
re-imaging optics were used, resulting in a plate scale of 
0.187" pixel^^ and a field of view of 3.'2. The seeing was 
~ 1". These observations consist of three 120 sec target- 
sky pairs, with the CCD dithered by 1.'5 between on-source 
positions. The NIR data are uncalibrated. 

The combined images were transformed to the World 
Coordinate System (and thereby to the same reference 
frame as the radio images) by registering to an image 
of the same area extracted from the Digitized Sky Sur- 
vey^ (DSS). The registration was accomplished by refer- 
encing the location of ^ 20 stars in common on both im- 
ages using the koords program in the Karma package 
(Gooch 1995) to perform a non-linear least squares fit for 
the transformation equations. The plate solution so found 
is accurate to a fraction of a pixel, and the overall regis- 
tration should be as good as that of the southern portion 
of the Guide Star Catalog, which is estimated to be 0'.'7 
(Taff et al. 1990). 

Deep optical images were constructed by applying a 9 x 9 
pixel boxcar median to the pixels with the lowest light 
levels, achieving limiting surface brightnesses of /iB=27.0 
mag arcsec"^ and ^r~2Q.?> mag arcsec"^. Color maps in 
B — R and B — K were made after convolving the im- 
ages to a common resolution, and using only pixels with 
a signal-to-noise greater than five in both bands. These 
maps will be shown here, but a full color analysis is de- 
ferred to a later paper (Evans et al., in prep.). Finally, 
a smoothed i?-band image was made by replacing stars 
with background values and convolving to the resolution 
of the intermediate resolution H i data. This image is used 
to evaluate the Hi mass to blue light ratio {Mhi/Lb) of 
various regions. 

3. RESULTS 

The large scale distribution of the H i is shown in Figure 
1^, where the left panel shows the low resolution H i con- 
toured upon the DSS image, while the right panel shows 
the intermediate resolution H i in greyscales with contours 
from the smoothed i3-band superimposed. Our new obser- 
vations delineate the distribution of tidal H i much more 
clearly than earlier H i observations, showing that the gas 
in the northern tail extends further than previously known, 
and revealing a small H i-rich disk companion to The An- 

^The Palomar Observatory Sky Survey was made by the California Institute of Technology with funds from the National Science Foundation, 
the National Geographic Society, the Sloan Foundation, the Samuel Oschin Foundation, and the Eastman Kodak Corporation. The Oschin 
Schmidt Telescope is operated by the California Institute of Technology and Palomar Observatory. 



tennae lying just beyond the southern tail. These results 
are confirmed by the ATCA observations by Gordon et al. 
(2001). Using the NASA Extragalactic Database (NED), 
this latter object is identified with the galaxy ESO 572- 
G045. We break the Hi emission into four separate com- 
ponents: the southern tail, the northern tail, the inner 
disks, and the new dwarf companion ESO 572-G045. The 
H I properties of each of these components are listed in 
Table | 

The large scale distribution and kinematics of the H I 
derived from the intermediate resolution data are shown 
in Figure I The upper left panel shows a false color rep- 
resentation of the data, with the H i shown in blue and 
starlight in green and white. The upper right panel of 
Fig. H shows the B — R color map with i?-band surface 
brightness contours superimposed. The lower left panel 
shows a color representation of the intensity weighted H I 
velocity field with iso- velocity contours superimposed. Fi- 
nally, the lower right panel shows a color representation of 
the H I velocity dispersion, again with H i column density 
contours superimposed. 

As revealed by the WSRT observations of van der Hulst 
(1979a), the majority of the Hi is associated with NGC 
4038 and the southern tail, with much less gas directly 
associated with NGC 4039 and the northern tail. We 
detect tidal H i emission over a velocity range of 1420 - 
1850 km s~^, similar to that found in the synthesis ob- 
servations of Mahoney et al. (1987), but less than found 
in the broader bandwidth observations of Gordon et al. 
(2001). We measure a total Hi flux of at least 54.5±0.5 
Jy km s~^ for NGC 4038/9, corresponding to a Hi mass 
of >4.7x10'^Mq (Table |). This is not very different from 
the total flux measurement of 57.1 Jy km s~^ measured 
at the ATCA by Gordon et al. (2001). However, when 
we divide the emission among the different components 
we do find differences. We measure a total flux of 49.7 
Jy km s^^ associated with the main disks and southern 
tail, compared to the 55.1 Jy km s~^ measured by Gor- 
don et al.. Examining the channel maps of Gordon et al., 
we find a mean flux level of ~ 10 mJy beam~^ (for a 40" 
beam) over the velocity range which we used for contin- 
uum subtraction. This suggest that, due to our improper 
continuum subtraction, we miss a total of 6 Jy km s^^ over 
the entire 600 km s~^ velocity range mapped by Gordon 
et al., which explains the different disk flux measurements. 
We also find differences for the northern tail and compan- 
ion galaxy, where we recover over twice the flux found by 
Gordon et al.. 

We next present details of the observations for each of 
the components individually, as well as for the tidal dwarf 
candidate(s) within the southern tail. 

3.1. Southern Tail Morphology 

The majority of the the tidal gas is associated with the 
southern tail. This tail extends to a projected radius of 
11^6 (65 kpc), and contains 2.8xl0^Mo of Hi along its 
entire length. The column density slowly decreases with 
increasing distance along the tail, reaches a minimum, and 
increases again near its end. Perpendicular to the tail, 
the gas is distributed more broadly than the optical light. 
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Fig. 1. — Left: Integrated Hi emission contoured upon on a greyscale representation of the DSS image of "The Antennae". The cross 
marks the phase center of the VLA observations, while the dashed circles represent the 75% and 50% response points of the VLA primary 
beam. The two main disks, NGC 4038 and NGC 4039, are labeled, as is the companion galaxy ESO 572 G-045. The contours are from 
the low resolution (40") datacube, and are are drawn at flux levels of (5, 10, 20, ...) X14.49 mjy beam~^ km s~^, with the lowest contour 
corresponding to an H I column density of 5x10^^ cm~^. Right: Greyscale representation of the intermediate resolution Hi (20.7" X 15.4"), 
with contours from the B-band image, after star subtraction and convolving to a resolution of 25". Optical contours are drawn from 27 mag 
arcsec"^ to 20.5 mag arcsec"-^ in intervals of 0.5 mag arcsec"^. The box-like optical contours to the northwest of the disk indicate the extent 
of the CCD frame in this direction. The scale bar indicates a length of 10 kpc for a distance of 19.2 Mpc. 



with a notable extension running along the outer (south- 
ern) edge of the tail. There is no optical counterpart to 
this outer extension at our limiting surface brightness of 
fiB < 27 mag arcsec"^ (see Fig. |^), with a resulting Hi 
gas-to-light ratio of Mhi/Lb> 1-5 MqLq^, compared to 
~0.5 Mq Lq^ on the optical tail. 

The gas within the southern tail shows considerable 
structure on the scale of one beamwidth (~ 1-2 kpc). This 
is further illustrated in Figure ^, which displays three or- 
thogonal projections of the intermediate resolution H i dat- 
acube. The upper left panel shows the sky view, I{a,6), 
and to either side of this are the two orthogonal position- 
velocity profiles, I{a,Vz) (bottom) and I{Vz,S) (right). 
The (a, S) view is displayed using a "maximum voxel" 
function. This function assigns to each pixel the maxi- 
mum intensity found along the velocity axis. The position- 
velocity plots are constructed in the more traditional man- 
ner, summing the emission along one of the spatial dimen- 
sions of the datacube, either right ascension or declina- 
tion. We will us e th ese plots below when examining the 
tail kinematics (§3.3) 



The maximum voxel display emphasizes regions where 
the H I is cold and dense (maximum number of H i atoms 
per unit velocity). As a result, it emphasizes the dynam- 
ically cold tidal tails. Note especially how the disk H i, 
which has the highest integrated H i column density, ap- 
pears fainter in Fig. ^ than in the zeroth moment maps 
(Fig. |l|). This is because the disk gas is spread over a 
very broad range of velocities. Fig. ^ reveals a number 
of dense knots within the southern tail, especially toward 
its end. The maximum voxel display ensures that this is 
not due to a line-of-sight integration effect where the tail 
bends back along our line-of-sight. Notice also that similar 
dense knots are not found within the northern tail. The 
southern tail also exhibits an interesting parallel or "bifur- 
cated" structure that starts where the southern tail begins 
to bend westward, and joins back together just before the 
location of the star forming regions identified by Mirabel 
et al. (see Fig. ^for the precise location of the star forming 
regions). 

The tail structure and kinematics are further illustrated 



in the channel maps, shown in Figure ^ (low resolution 
data after Hanning smoothing by a factor of two in veloc- 
ity to Au=10.6 km s~^, contoured upon the optical) and 
Figure || (intermediate resolution data; Av=5.2 km s~^, 
shown as a greyscale). Only channels containing emission 
from the tidal tails are shown in these figures, as the disk 
emission is spread over a much broader range of velocities 
and are shown separately below (§3^). Fig. Q gives a more 
complete mapping of the gas within the tidal tails and 
shows its relationship to the underlying starlight, while 
Fig. ^ allows a more detailed view of the structure of the 
gas within the tails. 

The greyscale representation of the channel maps in 
Fig. H clearly shows the "bifurcated" morphology of the 
southern tail mentioned above. The bifurcation is partic- 
ularly apparent in the channels between 1760 and 1734 km 
s^^. The two filaments merge together in a sideways "V" 
shape in the vicinity of the star forming regions identified 
by MDL92 (panel at 1734 km s^^ in Fig. |). The outer 
(southern) filament is associated with the previously de- 
scribed high Mhi/Lb material that lies off of the optical 
tail (see also panels at 1729 - 1772 km s^^ in Fig. |). The 
inner (northern) filament projects onto the optical tail and 
has a higher characteristic column density. 



3.2. Northern Tail Morphology 

The low resolution moment map (Fig. ^ shows that the 
northern H I tail extends beyond the end of the optical 
tail, to a projected radius of 7.'3 (40 kpc) from the center 
of NGC 4039. This extension is confirmed by the ATCA 
observations of Gordon et al. (2001). We detect a total of 
4.2x IO^A/q of H I associated with the northern tail, which 
is much more than found in any of the previous observa- 
tions. The northern tail is redder (Fig. 0b), less luminous, 
and has a lower relative H i content than the southern tail 
(Table |). 

We also notice that the gas in the northern tail is rele- 
gated to the outer half of the optical tail and in a purely 
gaseous extension beyond this. In particular, there is a 
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Fig. 2. — Large scale distribution of the H I and optical light in NGC 4038/9. The upper left panel shows a false color representation of the 
intermediate resolution H I data (blue) , and the CTIO B + V + R combined image in green and white. The upper right panel of shows the B — R 
color map, where colors have been determined for pixels with a S/N > 5 (p,B < 26.4 mag arcsec"^). Blue colors correspond to{B — R < 0.5], 
light blue/cyan to [0.5 < B - R < 0.7], green to [0.7 < B - R < 0.9], yellow/orange to [0.9 < B - R < 1.1], red to [1.1 < B - R < 1.3], 
magenta to [1.3 < B — R < 1.5] , and white to [B — R > 1.5] . Surface brightness contours from the star-subtracted B-band image convolved to 
25" resolution are drawn at intervals of (26.5, 26, 25.5, 20.5) mag arcsec"'^. The lower left panel shows the intensity weighted Hi velocities 
(green indicates velocities close to systemic, and red and blue indicate redshifted and blueshifted velocities, respectively) with isovelocity 
contours drawn at 10 km s~^ intervals. The lower right panel maps the Hi velocity dispersion with the intermediate Hi column density 
contours superimposed. In this panel blue/cyan corresponds to [aui < 5.0 km s~^], green to [5.0 < (Jhi < 9.0 km s~^], yellow/orange to 
[9.0 < afii < 13.0 km s~^], and red to [uhi > 13.0 km s~^]. Hi column densities are drawn at levels of (1, 2, 4, 8, ...) X28.88 mjy beam~^ 
km , where the lowest contour corresponds to 1x10^" cm~^. 
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II. 

Fig. 3. — Three orthogonal projections of the the intermediate resolution Hi datacube. The upper left panel shows the sky view, I{a, 5), 
and to either side of this are the two orthogonal position-velocity profiles, I{a, Vz) (bottom) and I{Vz, S) (right). The {a, 5) view is displayed 
using a "maximum voxel" function, which displays the value of the brightest pixel along the third (Vi) dimension (see text). The I{a, Vz) and 
I{Vz,&) projections are constructed by summing the emission along one of the spatial dimensions of the datacube, either right ascension or 
declination. These panels clearly show the well-behaved kinematic structure of the tidal features, as well as the wealth of gaseous sub-structure 
within the tails and disk. 



distinct gap in the H i distribution from a projected ra- 
dius of 25 kpc back to where the tail connects onto NGC 
4039. The difference in the H i content of the tails appears 
to be mirrored in the quantities of H i within the disks as- 
sociated with either tail (with NGC 4038 being gas rich 
and NGC 4039 being gas poor; see §3.5 below). These 
differences and the lack of H i at the base of the north- 
ern tail were already noted by van der Hulst (1979a), who 
attributed them to a difference in the H i content of the 



parent galaxies. We will return to this point in §4.1 



3.3. Global Tail Kinematics 

Velocities along the northern tail show a regular gradi- 
ent from 1610 km s~^ at the base of the Hi to 1540 km 
s~^ at the end of optical tail, while the velocities along 
the southern tail range from 1590, where the tail connects 
onto the northern disk, to a peak of 1770 km s~^ mid- 
way along the tail, declining to 1610 km s~^ toward the 
tip. This kinematic continuity is convincingly illustrated 
in the position-velocity plots in Fig. ^, which also demon- 
strates that the tidal gas is quite dynamically cold. The 
smooth run of velocities along the tails is already apparent 
from the Westerbork data of van der Hulst (1979a), and is 
expected from tidal interaction models (e.g. Barnes 1988, 
Hibbard & Mihos 1995). 

The well-behaved kinematics enable us to infer the ap- 
proximate tail geometry. For example, along the southern 
tail the H i column density is lowest at the regions of the 
most extreme redshift, and increases to either side of this. 
This indicates that this region of the tail lies perpendicular 
to our line-of-sight, with its motion directed most nearly 



toward us. At this location, the tail has a minimum pro- 
jected thickness and hence lowest H i column density. To 
either side of this the tail curves toward our line-of-sight, 
resulting in a larger pathlength through the gas-rich tidal 
material and hence a higher projected column density. The 
sharper velocity gradients to either side reflect the slew- 
ing of the velocity vectors as the tail bends away from us. 
Similarly, the sharp up-turn at the end of the southern tail 
in both space and velocity suggests that this region curves 
away from us. We will address the implications of this 



geometry for any embedded tidal dwarf galaxy in §3.4. 

Since tails are kinematically expanding structures 
(Toomre & Toomre 1972, Barnes 1988), the relative ve- 
locities of the two tails indicate that the southern tail is 
moving away from us and connects back to NGC 4038 from 
behind while the northern tail is swinging slightly toward 
us and connects back to NGC 4039 from the front. This 
geometry is supported by the lack of a dust absorption 
feature associated with the southern tail as it crosses from 
the south back to NGC 4038. In fact, the optical colors of 
these regions are quite blue (Fig. ^), whereas we would 
expect red colors from dust absorption if this gas-rich ma- 
terial lies between us and the disk. As shown directly by 
van der Hulst (1979a), these kinematics are consistent with 
the prograde spin geometry derived by the early numerical 
models of this system by Toomre & Toomre (1972; see also 
Barnes 1988; Mihos, Bothun & Richstone 1993). 

Following the southern tail kinematics back to its pro- 
genitor disk (NGC 4038, to the north), we see that the 
velocities become blueshifted toward the base of the tail 
(Fig. 0). Since the predominantly redshifted velocities of 
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Fig. 4. — H I channel maps of the tidal regions of NGC 4038/9 contoured upon the blue DSS image. The more sensitive low resolution Hi 
cube is used after smoothing by a factor of two in velocity (to a 10.6 km channel width). The 40" beam size is indicated in the lower left 
hand corner of the first panel, and each panel is labeled with its heliocentric velocity. Contours are drawn at levels of (3, 6, 12, 24, ...) xl.l 
mjy beam~i, where 1.1 mjy beam"^ is the single channel noise level, corresponding to a column density of 8x10^* cm~^. Only channels 
showing H l emission from the tails are shown. 
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Fig. 5. — A grayscale representation of the intermediate resolution H I channel maps of the tidal regions of NGC 4038/9. The greyscales 
are saturated at 1.1x10'^'' cm^^. Each panel is labeled with its heliocentric velocity. Only channels showing Hi emission from the tails are 
shown. This figure emphasizes the wealth of gaseous structure present in the tidal regions. 
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the tail constrain it to be swinging away from us and 
connect back to the disk of NGC 4038 from behind, the 
blueshifted velocities at the base indicate material that is 
streaming back toward the disk from the tidal tail (Hib- 
bard et al. 1994). Such streaming is a natural consequence 
of such encounters (e.g. Barnes 1988, Hernquist & Spergel 
1992, Hibbard & Mihos 1995), and arises because the tail 
material at smaller radii is more tightly bound than ma- 
terial further out. As a result, this material reaches its or- 
bital apocenter more quickly and subsequently falls back 
toward the inner regions. 

We can use the maximum projected tail length of 65 kpc 
and the disk rotation speed inferred from the Fabry-Perot 
observations of Amram et al. 1992 (155 km s~^, corrected 
for disk inclination), to infer an interaction age (measured 
from orbital periapse, when the tails were launched) of 420 
Myr. This is in excellent agreement with the best-match 
time of 450 Myr from the numerical model of NGC 4038/9 
run by Barnes (1988). It is also remarkably close to the age 
of the population of 500 Myr old globular clusters found 
by Whitmore et al. (1999; see also Fritze-v.Alvensleben 
1998), supporting the hypothesis that this GC population 
formed around the time when the tails were first ejected 
(Whitmore et al. 1999). 

From Fig. ^ we note that the velocity dispersion within 
the tidal tails is actually quite low: except for the base of 
the southern tail, the H i dispersion is less than 16 km s^^, 
and the mean of 7 km is typical of values measured 
in undisturbed disk galaxies. This is true of long-tailed 
mergers in general (Hibbard & van Gorkom 1996, Hib- 
bard & Yun 1999a). The dispersion near the base of the 
southern tail has higher characteristic values of 13-25 km 
s~^. These regions likely lie along a steeper gradient in 
the potential, and are also regions where the line-of-sight 
velocities are changing rapidly (Fig. ||:). These effects will 
lead to larger real and apparent velocity gradients, increas- 
ing the single-beam dispersion. 

3.4. The Tidal Dwarf Galaxy Candidates 

The improved resolution of the VLA observations al- 
lows us to directly address the behavior of the atomic gas 
in the vicinity of the putative tidal dwarf galaxy(s) iden- 
tified near the end of the southern tail of The Antennae. 
We show a close-up of this area in Figurco, which has the 
same arrangement as the panels in Fig. T 

Schweizer (1978) was the first to identify four Hii re- 
gions and blue U BV colors associated with the end of the 
southern tail. He also identified a patch of low surface 
brightness material bending sharply to the north after the 
end of the optical tail, which he suggested was a sepa- 
rate dwarf stellar system. We indicate the approximate 
location of this candidate TDG in Fig. ^ using a thick 
yellow contour drawn at approximately one half of the lo- 
cal peak H i column density after background subtraction 
and labeled TDG [S78]. Schweizer noted that in the low- 
resolution Hi map of van der Hulst (1979), the Hi ap- 
pears to be more closely concentrated in the low surface 
brightness extension of the tail, and hypothesized that the 
dwarf may have been created during the interaction, as 
envisioned by Zwicky (1956). 

The end of the tail was subsequently studied by Mirabel 
et al. (1992), who concentrated on what they called a "de- 



tached condensation of gas and stars at the tip of a tidal 
tail", consisting of a twisted stellar bar embedded in an 
envelope of diffuse optical emission and containing three 
emission line complexes. We indicate the approximate lo- 
cation of this candidate TDG in Fig. || using a thick white 
contour drawn at approximately one half of the local peak 
H I column density after background subtraction and la- 
beled TDG [MDL92]. The three star forming regions iden- 
tified by MDL92 are labeled as Regions I— III in Fig. ||a, 
and are represented by black or white dots in the remain- 
ing panels. 

While all subsequent observers have referred to the re- 
gion identified by MDL92 when discussing the putative 
TDG, we will address the dynamical nature of both of 
these regions. This is not a trivial task, as the precise 
boundaries of the TDG candidates have never been ex- 
plicitly defined. 

Our H I observations (as well as those by Mahoney et al. 
1987 and Gordon et al. 2001) are a vast improvement on 
the observations of van der Hulst (1979), and show that 
the tail is a continuous kinematic structure. Specifically, 
there is no "detached" region which can readily be identi- 
fied as a distinct entity. Further, the velocity field varies 
smoothly through the region containing the high gas col- 
umn densities and star forming regions without a signif- 
icant twist or kink. The strongest velocity gradient cor- 
responds with the region where the tail turns abruptly to 
the north (compressed contours near the center of the map 
in Fig. with a similar gradient just to the east of this). 
This smooth gradient suggests that the velocity field in 
this region is dominated by projection effects, with the 
tail bending away from us back into the plane of the sky; 
as a result the velocity vector slews from pointing toward 
us to pointing away from us. Therefore, we find no clear 
kinematic evidence in the velocity field for any dynami- 
cally distinct entities. It is possible that the strong tidal 
gradient may mask the kinematic signature of any mass 
concentration. 

Since there are no well defined boundaries to the dwarf 
candidates, we calculate the H i mass and blue luminosity 
in their vicinity by summing the emission within succes- 
sive circular apertures centered on the local H i column 
density peaks. The maximum aperture radius is set by the 
half-width of the tail. A mean background is subtracted, 
and the H i flux is summed within each aperture of radius 
R. The optical luminosity is measured from the smoothed 
star-subtracted i?-band within the same apertures. 

The results of these calculations are shown in Fig. ^ 
The left panel shows the result for TDG [MDL92] and the 
right panel shows the result for TDG [S78]. In these plots, 
the solid and dashed lines show how the H i mass and blue 
luminosity grow as a function of the aperture radius R. In 
Table |^ we list the total H i mass and optical luminosity 
measured by these curves for each dwarf candidate, but 
since there are no distinct boundaries for either region, 
these values should be considered very rough guides. 

The strong velocity gradients mentioned above give rise 
to the regions with the highest H i velocity dispersions in 
Fig. as gas with different space velocities fall within 
the same beam. Aside from these regions, there are other 
local peaks in the gas velocity dispersion which might pro- 
vide possible evidence for a mass concentration: between 
Regions I & II and Region III the Hi velocity disper- 
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Fig. 6. — Four-panel mosaic of the region of the putative Tidal Dwarf Galaxy. The upper left panel shows the optical image in false 
color, with contours from the intermediate resolution Hi data drawn at levels of 14.44 Jy km s~^ (5x10^^ cm~^) xl.5", n = 1 — 9. These 
contours are repeated in all panels. The three star forming regions identified by Mirabel et al. (1992) are labeled in the upper left panel, and 
are indicated by black or white dots in subsequent panels. Solid heavy contours delineate the approximate half-light level around the Tidal 
Dwarf candidates identified by Mirabel et al. 1992 (heavy white contour labeled TDG [MDL92]) and by Schweizer 1978 (heavy yellow contcuir 
labeled TDG [S78]). The upper right panel shows a color representation of the B — R color map, with the same color mapping as in Fig. Hb. 
The lower left panel shows the intensity weighted Hi velocities with isovelocity contours drawn in black at 5 km s~^ intervals. The lower 
right panel maps the H I velocity dispersion, with the same color mapping as in Fig. Hd. 




Radius (kpc) Radius (l^pc) 



Fig. 7. — Results of luminous mass calculation for the candidate Tidal Dwarfs identified by Schweizer 1978 (right panels, labeled TDG 
[S78]) and by Mirabel et al. 1992 (left pappls, labeled TDG [MDL92]). On both cases, the origin is taken to be the peak gas surface density 
in the integrated intensity map (see Fig. p|) . The two curves show the enclosed H I mass (solid curve) and optical luminosity (dashed curve) 
as a function of the circular aperture radius R. 



sion increases to ^12 km s~^ vs. an average of 6-7 km 
s~^ along the tail. Similar signatures are seen at the lo- 
cation of tidal dwarf candidates within the optical tails of 
the merger remnants NGC 7252 (Hibbard et al. 1994) and 
NGC 3921 (Hibbard & van Gorkom 1993, 1996). How- 
ever, it is noteworthy that the three H ii regions actually 
fall on the edges of H i maxima. Similar signatures are 
seen near giant H il regions in dwarf galaxies (e.g. Stewart 
et al. 2000, Walter & Brinks 1999), where the increased H i 
linewidth is due instead to kinetic agitation of the gas from 
energy deposited by young stars, SNe and stellar winds 
(e.g. Stewart et al. 2000, Yang et al. 1996, Tenorio-Tagle 
& Bodenheimer 1988 and references therein) rather than 
the gravitational effects of a mass concentration. 

There is also an increased dispersion to the north of 
TDG [S78] (dispersion increases to 15.6 km s^^ over the 



surrounding value of 11 km s~^). However, an examina- 
tion of the line profiles shows that the increased dispersion 
comes from low-level emission spread over many channels. 
Adjacent regions have this broad low- level component, but 
also have a brighter narrower component which leads to a 
lower intensity-weighted line-width. So we do not believe 
that this is a signature of a mass concentration at this 
location. 

Since the moment maps are susceptible to line-of-sight 
integration effects, we turn to the velocity cubes directly 
for further insight. From the channel maps, we find that 
the gas near the TDG obtains the highest gaseous phase- 
space densities in this system (i.e. density enhancements 
localized in space and velocity). In fact, this region has 
a higher H i column density per unit km than any 
other region in this system, including regions within the 
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disk or overlap region (see also Fig. ||a). The peak of 
21 niJy beam^^ in a single velocity channel of width 5.2 
km s~^ corresponds to a single channel column density 
of 3.8x lO^*' cm~^, and the integrated column density of 
1.3xl0^^cm~^ represents an enhancement of 6 over the 
mean tail density. The high-resolution data shows an even 
higher single-channel peak of 6.4x10^° cm~^. 

A closer examination shows that there are actually two 
dense knots in this region. These are seen in the channel 
maps from the intermediate resolution datacube, plotted 
at full velocity resolution in Fig. H. The densest concen- 
tration appears at V=1708 km s and lies just north of 
star forming Region III. The second concentration en- 
compasses star forming Regions I & II and appears at 
y=1719-1729 km s~^. Such concentrations suggest re- 
gions of gaseous dissipation. While this is not unequivocal 
evidence of a self-gravitating dwarf-sized object, since in- 
creased gaseous dissipation is expected in regions of recent 
star formation, it does indicate that there is something 
unique about this region of the tail. 

The channe l m aps also show the bifurcated structure 
mentioned in § ^.l| . It is very interesting that the two par- 
allel filaments join just east of the location of the star 
forming regions associated with the TDG (the "V" -shaped 
feature seen in the panel at 1734 km s^^ of Fig. ^). Again, 
these observations suggest that there is something special 
about the gas and tail geometry at the location of the star 
forming regions. 

The parallel filaments can be seen separately in the 
channels from V= 1755- 1724 km s"^ The northern fil- 
ament quite clearly lies along the optical tail, while the 
southern filament is displaced by about 8 kpc to the south. 
The density peaks are of similar magnitude in either fil- 
ament, but only at the location of the putative TDG is 
there a corresponding increase seen in the optical light 
underlying the gas peaks. The gas concentration associ- 
ated with Region III appears to lie within the northern 
filament, while it is difficult to say whether the concentra- 
tion associated with Regions I & II are associated with 
one filament or the other. At velocities below 1693 km 
s~^, it appears that only the southern filament continues 
on to the gas rich "hook" at the end of the tail. 

Fig. 1^ shows that the blue tail colors discovered by S78 
are not restricted just to the high column density gas near 
either TDG, but represents a trend of increasing blueness 
with distance along the tail. We find that while the outer 
disks and southern tail have similar B — R colors, the end 
of the southern tail is bluer \n.V — R than B — V while 
the reverse holds for the outer disks. This suggests that 
the tail material is younger and/or more metal poor than 
the outer disk material (Schombert et al. 1990, Weilbacher 
et al. 2000), although we postpone to a later paper a full 
color analysis (Evans et al. 2001). 

3.5. Disk Morphology and Kinematics 

Fig. H shows an overlay of the high resolution H i con- 
tours upon a false-color representation of the i?-band im- 
age (Fig. ^), of the B - K color map (Fig. ||b), of the H i 
velocity dispersion map (Fig- ^), and of the NIR K-ha.nA 
image (Fig. ^). In Fig. ^ we show a color representation 
of the H I velocity field with isovelocity contours drawn 
at intervals of 5 km s^^. Finally, in Fig. ^f we show the 



integrated CO map of Wilson et al. (2000) contoured in 
green upon a greyscale representation of the H i column 
density. The latter image gives a better idea of the local 
minima and maxima in the H i distribution, as these are 
not always self-evident in the contour map. In this panel 
red crosses indicate the location of the "Super Star Clus- 
ters" (SSCs) identified in the HST image of Whitmore et 
al. (1995, from their Table 1). 

In Fig. ^b dark regions represent the reddest B ~ K col- 
ors, generally indicating regions of high dust extinction, 
while light regions represent the bluest colors, generally 
indicating unobscured massive young stars. Comparison 
of Figs. 0b & f shows that the CO is closely associated 
with the red features in the B — K map, so that the CO 
traces the dust distribution. 

Figs. ^ & f show that on a global scale, the H i distribu- 
tion traces both the regions containing the most star clus- 
ters and the regions with the highest dust content, e.g. the 
star forming ring in NGC 4038 and the star forming half 
of the disk of NGC 4039 (to the northeast of the nucleus), 
and the dusty region between the two galaxies and to the 
southwest of NGC 4039. On finer scales (i.e. of order a few 
kpc) , the H I is not particularly associated with the bluest 
star clusters or reddest dust concentrations. Most of the 
SSCs in Fig. ^ actually appear to lie along steep gradients 
in the H i and CO distributions. The molecular gas peaks 
are often displaced from the Nhi peaks and many of the 
SSCs actually lie in the transition zone between the two 
gas phases. Where there are no displaced CO peaks, the 
SSCs still lie predominantly along steep gradients in the 
H I distribution. 

The NIR image in Figs. ^ does not show additional 
light peaks in the regions of highest Hi column density 
(recall that dust extinction is 10 times weaker in the NIR 
compared to the optical), so we conclude that the occur- 
rence of SSCs along gradients in the H i distribution is not 
due to dust associated with the H i obscuring clusters at 
the regions of highest H i column density. Rather, we sus- 
pect that it is due to young stars ionizing or sweeping the 
H I on the one hand, and the conversion of H i to molecular 
form on the other. We refer the reader to Zhang, Fall & 
Whitmore (2001) for a much more thorough discussion of 
the relationship between SSCs of various ages and the H i. 
Radio continuum, CO, X-ray, FIR, MIR, Ha and optical 
morphology and the H i and CO kinematics. 

The lower H i contours in Fig. ^ are very close to- 
gether, with the H i column density falling off sharply be- 
low 4x10^° cm~^. The western part of the disk of NGC 

4038 is largely devoid of H i, as is most of the extended 
disk of NGC 4039. Examination of the lower resolution H i 
maps (Figs. |[ ||a) shows that this is true even at lower H I 
column densities. Probably not coincidently, these regions 
lack current star forming regions. 

The intensity weighted Hi velocities of NGC 4038/9 
(Fig. He) shows a very disturbed velocity field, similar 
to those mapped in Hi (Gordon et al. 2001) Ha (Bur- 
bidge & Burbidge 1966, Rubin et al. 1970, Amram et al. 
1992) and in CO (Stanford et al. 1990, Wilson et al. 2000, 
Zhu 2001). In particular, the disk of NGC 4038 has a 
north-south gradient spanning ~ 250 km s^^, while NGC 

4039 has a northeast-southwest gradient spanning a sim- 
ilar range (see also Fig. The general sense of motion 
— north blueshifted, south redshifted — agrees with the 
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Fig. 8. — Hi channel maps of the region around the putative Tidal Dwarf Galaxy within the southern tail of NGC 4038/9 contoured upon 
a greyscale representation of the B-band image. The greyscales range from /is=28.5 mag axcsec"^ (white) to 25 mag arcsec"^ (black). The 
intermediate resolution Hi cube is used at its full velocity resolution (5.2 km channel width). Contours are drawn at levels of 2.8 mjy 
beam~^ x 1.5", n = — 4, where 2.8 mJy beam~^ corresponds to a single channel column density of 5x10^® cm~^. 
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Fig. 9. — Details of the Hi emission within the disks of NGC 4038/9. The upper left panel shows contours from the high resolution Hi 
datacube (11.4" X 7.4") upon a false-color representation of the B-band image. Contours are drawn at levels of (1, 2, 4, 8, ...) x7.74 mjy 
beam~^ km s~^, where the lowest contour corresponds to a column density of IxlO^'^cm"^. The upper right panel shows the same Hi 
contours upon the B — K color map, where lighter shades indicate bluer colors and darker shades represent redder colors. The middle two 
panels show the intensity weighted Hi velocities with isovelocity contours drawn at 5 km s~^ intervals (left) and the Hi velocity dispersion 
with the same Hi contours as in the first panel (right). The next two panels shows K-hand image with the Hi contours superimposed on 
the left, and the a greyscale representation of the Hi image on the right, with the CO contours of Wilson et al. (2001) superimposed and red 
crosses marking the location of the Super Star Clusters identified by Whitmore & Schweizor (their Table 4). 




sense of rotation inferred from the tidal tail. However, be- 
cause of the extremely disturbed nature of both disks, the 
kinematics of these regions are not simply rotational. It 
is likely that there are strong streaming motions along the 
star forming loop and along the overlap region connecting 
the two disks. 

In Figure |l^ we present the disk channel maps, made 
from the intermediate resolution data after Hanning 
smoothing by a factor of four in velocity (Au=20.4 km 
s""'^). Instead of a normal rotational pattern, the redshifted 
gas concentration to the SW of NGC 4039 shows very lit- 
tle spatial gradient over the velocity range V=1836-1708 
km s~^, while the blueshifted gas to the northeast of NGC 
4039 has very little spatial gradient over the velocity range 
V=1579-1472 km s~^. The gas not associated with ei- 
ther of these two regions appears to flow from the NE at 
V=1558 km s"! to the SW at V=1665 km s'K This is 
likely material which is flowing along the overlap region 
between the two galaxies. The H i in the northeast part of 
NGC 4039 may be material which is being accreted from 
NGC 4038 along the bridge connecting the two systems. 



We also notice that the disk gas connects smoothly to the 
gas associated with the base of the southern tail (chan- 
nels at V=1472-1686 km s~^). From the tail kinematics 



(§|3.3D we concluded that the blueshifted velocities associ- 
ated with the base of the tail require that this gas streams 
back onto the outer disk of NGC 4038 from the tail. In 
Fig. ^ we find a region of high velocity dispersion along 
the northern loop of gas (dispersion of 28 km s~^ vs 22 km 
s~^ at adjacent regions). We suggest that this increase is 
associated with gas falling back from the tail mixing along 
the line-of-sight with gas associated with the disk. 

3.6. The Nearby Dwarf Galaxy ESQ 572-G045 

As part of this observation we also detect Hi associ- 
ated with the Scd dwarf galaxy ESO 572-G045 just to the 
southwest of the tip of the southern tail (5.7' or 32 kpc 
from the TDG). This system lies at a projected distance of 
90 kpc from The Antennae. This galaxy is very flattened, 
with an axis ratio near 10, associating it with the class of 
"superthin" galaxies (Goad & Roberts 1981), and qualify- 
ing it for membership in the "Flat Galaxy Catalogue" of 
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Fig. 10. — Hi channel maps of the disk regions of NGC 4038/9 contoured upon the CTIO B-band image. The intermediate resolution 
Hi cube is used after smoothing by a factor of four in velocity (to a 21 km s^^ chaimel width). Each panel is labeled with its heliocentric 
velocity. Contours arc drawn at levels of (-3, 3, 6, 12, 24, ...) xO.4 m.Jy beam^^, where 0.4 mJy bcam^^ is the rebinned channel noise level, 
corresponding to a column density of 3x10^'^ cm^"^. These maps show the entire velocity spread of the present observations, and the reader 
should refer to Gordon et al. (2001) for H i emission from the disk regions that falls outside of our bandpass. 
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Karachentsev et al. (1993). It is curious that a similar su- 
perthin galaxy was also serendipitously discovered by Due 
et al. (2000) in their HI mapping of the interacting system 
NGC 2992/3, another interacting pair with a candidate 
Tidal Dwarf within a gas-rich tidal tail. 

Contours of the H i column density from the various res- 
olution data cubes are shown on the optical image from 
the DSS in Figure |ll|. We detect a total 3.7xlO*Mo of 
H I spread over 100 km s~^ and extending to a radius of 6 
kpc (i.e., twice as far as the optical light). The galaxy is 
nearly edge-on and shows a linear velocity gradient. As- 
suming an edge-on orientation and circular velocities (with 
Vcirc = 50 km s~^), this suggests a dynamical mass {Mdyn) 
of S.SxIO^Mq out to 6 kpc. Using the optical properties 
given in NED (major and minor axis of 1.01' x 0.11' and 
nriB =16.70 mag), we derive dimensions of 5.9 kpc x 0.6 
kpc and a luminosity oi Lb = 1.9 x 10^ Lq. This sys- 
tem is quite gas rich {Mni/Mdyn — 0.11, Mhi/Lb — 1.9), 
with global properties similar to those of other superthin 
galaxies (Matthews et al. 1999). 

The low resolution H i map in Fig. 



11a, clearly shows 



some low column density gas extending to the east. Since 
this extension is one-sided in nature, it is more indicative 
of ram pressure effects than tidal effects. The very thin 
and flat disk of ESO 572-G045 also argues against signif- 
icant tidal forces acting on it. This suggests that there 
may be an extended warm or hot gaseous halo around 
NGC 4038/9 that is affecting the diffuse cold ISM of ESO 
572-G045 (e.g. Moore & Davies 1994). 

4. DISCUSSION 

The wealth of information provided by these observa- 
tions is a testament to the utility of H i spectral line map- 
ping of disturbed systems. While in some sense the present 
data are just the beginning point of many more involved 
investigations, particularly for dynamical modeling of this 
system, they also lead to many interesting results on their 
own, which we discuss in this section. 

4.1. The lack of Hi at the base of the Northern Tail 

The difference between the northern and southern tail in 
both optical structure and H i content was already noted 
by van der Hulst (1979a), who attributed it to a difference 
in the H i content of the parent galaxies. The lack of H i 
within the main disk NGC 4039 (i.e. the progenitor of 
the northern tail) was seen as support for this interpreta- 
tion. However, we now feel that the observations are at 
odds with this explanation. For while it is true that the 
outermost regions of tidal tails should originate from the 
outermost regions of the progenitor disks, it is not true 
that the outermost regions of the progenitor disks end up 
only at the outermost regions of tidal tails. In fact, this 
gas rich material should extend all the way back to the 
host system. This fact is apparent already in the simple 
models and illustrations of Toomre & Toomre 1972 (see 
their Fig. 2). This reflects the fact that tails are not linear 
structures, but actually two dimensional "ribbons" twist- 
ing through space, with the outer disk material forming 
a "sheath" around the inner disk material (Mihos 2001). 
As a result, the gas rich outer regions should extend along 
the entire length of the tidal feature back toward the pro- 
genitor disk. It therefore seems that a low H i gas content 
alone cannot explain the gas-rich outer tail and gas-poor 



inner tail. 

In a recent paper, Mihos (2001) describes a kinemati- 
cally decoupling between the gaseous and stellar compo- 
nents of a tail that occur during tail formation. How- 
ever, this decoupling leads to a displacement between the 
gaseous and stellar tidal features, and does not totally re- 
move large quantities of gas from along radial segments 
of the tail. It is unlikely, therefore, that this mechanism 
explains the present observations. 

Hibbard, Vacca & Yun (2000) investigate several pos- 
sible mechanisms for creating differences between H i and 
optical tidal morphologies. Two plausible mechanisms are 
ram-pressure stripping of the tidal gas from an expanding 
superwind and ionization of the gas by UV photons escap- 
ing from the starburst. While there is some evidence for a 
nuclear outflow in NGC 4038/9 (Read, Ponman, & Wols- 
tencroft 1995, Sansom et al. 1996, Fabbiano, Schweizer & 
Mackie 1997, Lipari et al. 2001, Fabbiano, Zezas & Mur- 
ray 2001), the soft X-ray morphology (Read et al. 1995, 
Fabbiano et al. 2001) does not suggest that this wind is di- 
rected toward the northern tail, so we concentrate instead 
on possible ionization effects. 

The situation explored by Hibbard et al. (2000) is that 
the starburst has cleared enough dust and gas from the in- 
ner regions to provide direct unobscured sightlines to the 
tails. By equating the surface ionization rate with the re- 
combination rate, they find that the gas should be ionized 
out to a distance given by: 



< 25 kpc 



fe. 



1/2 



0.05 

2 X 102"cm-2 



1/2 



HI 



dL 
5 kpc 



1/2 



(1) 



Where Nhi and dL are the column density and thick- 
ness of the H I layer being ionized, fesc is the fraction of 
ionized photons emerging from the starburst, and we have 
made the standard assumption that most of the starburst 
luminosity is emitted in the far infrared (see Hibbard et al. 
2000 for details). 

The fiducial values used in eqn. |^ are appropriate for 
NGC 4038/9, and we thus find that the starburst may 
well be capable of ionizing gas out to the H i base of the 
northern tail (25 kpc projected radius), providing the tail 
has a clear sightline to the ionizing stars and that fesc is 
not much smaller than 0.05. Evidence that fesc is at least 
as large as this is provided by numerical work by Dove, 
ShuU & Ferrara (2000), who derive fesc ^ 0.1 for normal 
galaxies, and suggest that starbursts should have signifi- 
cantly higher values. Hurwitz, Jelinsky, & Dixon (1997) 
find upper limits to the escape fraction of UV photons 
of 3-57% for four IR luminous starbursts, while Bland- 
Hawthorn & Maloney (1999) calculate fesc = 0.06 for the 
Milky Way (see discussion in Bland-Hawthorn & Putman 
2001). 

So the question becomes, is it likely that the northern 
tail of NGC 4039 has a relatively clear view of the star- 
forming disk, while the southern tail does not? We believe 
the answer to this question is yes. The majority of the 
unobscured massive star forming regions in NGC 4038/9 
are associated with the disk of NGC 4038 and the mate- 
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Fig. 11. — Contours of the integrated Hi emission of the small superthin galaxy ESO 572-G045 overlaid on the DSS image. These panels 
show the Hi data from the low— (left), intermediate— (middle) and high- (right) resolution datacubes, with the beam size of each dataset 
drawn in the lower left corner. Contours are drawn at levels of [144.9 Jy km s"! x2"/2, 57.75 Jy km s"! x2"/2_ ig.i jy g-l x2"] 
(n = 0,1,2...), for the low, intermediate and high-resolution data, respectively. The lowest contour corresponds to a column density of 
[1 X lO^'' cm~^, 1x10^" cm~^, 4x10^*^ cm~^], respectively. Notice the gaseous extensions extending to the cast in the intermediate and low 
resolution maps. 



rial along the bridge (Whitmore & Schwcizer 1995, Whit- 
more et al. 1999). Since tidal tails are spun off close to the 
spin plane of the progenitor (Toomre & Toomre 1972), the 
southern tail will not have a clear view to the star forming 
regions located within the disk NGC 4038, and the body 
of NGC 4039 would block ionization originating from the 
bridge star forming regions. The body of NGC 4039, on 
the other hand, appears to have a perfect view of the star 
forming disk of NGC 4038 and the northern tail should 
have a view to the backside of this disk. This prediction 
can be tested by numerical models to constrain the space 
geometry of the tail. 

We suggest that NGC 4039 had a normal H i distribu- 
tion, but that much of this gas has been photoionized by 
the star burst. We believe a similar effect is responsible 
for the lack of H i at smaller radii in the otherwise gas-rich 
tidal tails of the NGC 7252, Arp 105, and Arp 299 systems 
(Hibbard, Vacca & Yun 2000). If this is indeed the case, 
then the ionized gas should be visible via its recombina- 
tion radiation. Using the equation given in Hibbard et al. 
(2000), the expected emission measure of this radiation 
should be EM — 0.8cm^^pc for a gas column density of 
2xl0^''cm^^ and thickness of 5 kpc. This is within the 
capabilities of modern CCD detectors (e.g. Donahue et al. 
1995; Hoopes, Walterbos & Rand 1999). 

4.2. The Bifurcation of the Southern Tail 

As mentioned above, tails are two dimensional "rib- 
bons" twisting through space. It is possible that a lat- 
eral twist may cause the outer edge of the ribbon to lie in 
a different plane from the inner regions, i.e. for the most 
gas- rich (highest Mhi/Lb) regions to appear in a different 
plane from the less gas-rich but optically brighter regions 
(Mihos 2001). In Hibbard & Yun (1999a) we suggest that 
this this effect may be exacerbated by a pre-existing warp 
in the progenitor disk, providing even more marked off- 
sets. This provides a simple explanation for the gas-rich 
but optically faint outer edge of the southern tail. Similar 
bifurcated tails are seen in several systems (e.g. M81, van 
der Hulst 1979b, Yun et al. 1994; NGC 3921, Hibbard & 
van Gorkom 1996; NGC 2535/6, Kaufman et al. 1997; Arp 
299 Hibbard & Yun 1999a). 

What is intriguing is that the two filaments appear to 
join back together just at the location of the star forming 
regions associated with the candidate tidal dwarf TDG 
[MDL92] (the "V"-shape seen in Figs. | & | and discussed 
in §0 & §3.4). We do not know whether this is a coinci- 



dence or an important clue. Further insight into this ques- 
tion will have to await a more detailed numerical study of 
this system. 

4.3. The Formation of the Tidal Dwarf Galaxy 

We wish to take advantage of the higher spatial and 
velocity resolution of our H i observations to examine the 
dynamical nature of material in the vicinity of the TDG 
candidates. This is a difficult task, since there are no dis- 
tinct objects which can be defined unambiguously. Fur- 
ther, if there are objects within the tail, they are clearly 
not in isolation, and the kinematics may not be due to 
the local mass concentration. All of these factors make a 
clear dynamical analysis fraught with uncertainties. Still, 
the present data are the best available, in terms of both 
spatial and kinematic resolution, and this is the nearest 
system with a candidate tidal dwarf, so we will proceed 
with all due caution. 

The first hypothesis we will test is whether there are 
concentrations of gas and light centered on the gas peaks 
with enough mass in luminous matter alone to account 
for the measured H i line-width. To this end, we compare 
the virial mass inferred from the Hi linewidth {Myir) to 
the total luminous mass inferred from the H i flux and the 
optical luminosity (M/„m). A value of Mium/Mmr — 0.5 
is sufficient for a region to be bound, while a value of 
Mium/Myir = 1 is required for the region to be virialized. 

To calculate the virial mass, we must make some as- 
sumptions about the 3-dimensional structure of any mass 
concentration. For the sake of simplicity, we assume that 
there is an isotropic spherical mass concentration of con- 
stant density. For this situation, the virial mass is given 
by Alvir — 3(7^ X arh / G, where G is the gravitational con- 
stant, a is the one-dimensional velocity dispersion, is 
the projected half-light radius, and a is a geometric factor 
(Binney & Tremaine 1987) ~ 2.74 for the adopted geom- 
etry. We take the velocity dispersion at the location of 
the peak column density, and use the deconvolved half- 
light radius found from fitting a Gaussian plus constant 
background to the H i mass distribution. 

The luminous mass is calculated from the observed mass 
and optical luminosity following the prescription given in 
Paper II. In particular, the possible contribution to the lu- 
minous mass by stars is estimated by adopting a range of 
stellar mass-to-light ratios (M*/Lb) representing different 
star formation histories (SFHs). The gas mass is derived 
from the Hi mass by multiplying by a factor of 1.36 to 



18 



H I in The Antennae 



E 0.6 



0.2 



: MOLQ^' ■ ' ■ 1 ■ ■ ■ ■ 


' ' ' ' . 




1 1 1 1 1 1 1 1 1 1 r 1 1 1 1 1 1 1 1 

: S78 : 







1 

t 0.8 
E 0.6 


^ 


" ■ t./L,-i'. 






a" 

— 0.4 






" 




u 

o 

0.2 


■ 


■~— <.A,=0; 

" 



23456 23456 

Deconvolved Radius (kpc) Deconvolved Radius (kpc) 

Fig. 12. — Results of dynamical analysis for the candidate Tidal Dwarfs identified by Schweizer 1978 (right panels, labeled TDG [S78]) 
and by Mirabel et al. 1992 (left panels, labeled TDG [MDL92]). On both cases, the origin is taken to be the peak gas surface density in the 
integrated intensity map (see Fig. BI). We plot Mi^^/M^yn as a function of the deconvolved aperture radius r. Three curves are drawn, each 
one representing a different stellar mass-to-light ratio: the thick solid curve shows the result for Mt/Lg = Mq Lq^, the dotted curve shows 
the results for M^/Lb = 2 Mq L'^ , and the dashed line shows the results for M»/Ls = 5 Mq Lq^. 



take into account a primordial abundance of Helium. No 
correction is made for the presence of molecular gas since 
recent CO observations at the location of the TDG by 
Gao et al. (2001) and Braine et al. (2001) find only trace 
amounts of molecular gas in the vicinity of the tidal dwarf 

We evaluate three different stellar mass-to-light ratios, 
which span a range of star formation histories. We eval- 



uate M^/Lb = MqL~^ (i.e., Ml, 



also ap- 



propriate for a young instantaneous burst); M^^/Lb = 2 
MqLq^, corresponding to a constant SFH at an age of 

10 Gyr, and M^/Lb ~ 5 MqLq^, representing an expo- 
nentially decreasing SFH with a time constant of 4 Gyr 
at an age of 10 Gyr (see Paper II for more details). The 
results of the luminous to virial mass calculation are given 
in Table |. 

The ratios of the luminous to virial mass given in Ta- 
ble ^ are less than one, suggesting that neither mass con- 
centration is virialized. However, they do not rule out the 
presence of a bound object. They do imply that there is 
not a sufficient gas density for the regions to be bound 
by gas alone (i.e., Mium/Mmr < 0.5 for M^/Lb — 
Mq Lq^). This contradicts the conclusion reached for this 
same region by Braine et al. (2001). The reason for this 
discrepancy is that Braine et al. use the much narrow CO 
linewidth to calculate a dynamical mass, yet include both 
H I and stellar light in their luminous mass calculation 
(as mentioned above, the mass in molecular mass inferred 
from the CO line flux is negligible compared to the Hi). 
The fact that the H i linewidth is much broader than the 
CO linewidth argues that a much larger mass concentra- 
tion is required to bind it, and the Braine analysis is inap- 
propriate. This point is discussed in more detail in Paper 
H. 

A further conclusion from the above analysis is that if 
the dwarf candidates are bound, they require a signifi- 
cant contribution to the total mass from evolved stars, or 
they are bound by dark matter. In this case, the calcu- 
lated values of Mium/Mmr are valid only if the H i line- 
width also characterizes the motion of the stars (and/or 
dark matter). This is probably not the case, since in nor- 
mal disks the older stellar populations (i.e., those with 
higher M*/Lb ratios) have a signiflcantly higher disper- 
sion than the gas, by factors of 2-3 (van der Kruit 1988), 
which would lead to a virial mass estimate a factor of 4- 



10 higher and Mium/Mmr ratios correspondingly smaller. 
Young stellar populations have a velocity dispersions sim- 
ilar to the cold gas, but also have smaller M^/Lb ratios. 
A measurement of the stellar kinematics and a better es- 
timate of the stellar mass-to-light ratio is needed to be 
more conclusive. Finally, we note that the expected cir- 
cular velocity given by = GMium/R is 20-30 km s^^. 
Fig. ^ shows no sign of rotation at this level. We con- 
clude that either TDG candidate may be, but need not 
be, self-gravitating. 

This virial analysis includes many caveats. Primar- 
ily, the above analysis implies that we know the ratio 
Miurn/Myir to within a factor of two, and in truth we do 
not know either quantity to this precision. Further, unless 
the bound regions are much smaller than the width of the 
tail, the geometry is unlikely to be spherical, and the ve- 
locities are unlikely to be isotropic. Finally, any entrained 
objects are clearly not in isolation, complicating the virial 
analysis. We will be addressing these uncertainties using 
simulated observations of bound objects within N-body 
tails (Kohring, Hibbard & Barnes, in preparation). 

The above analysis depends on a specific geometry and 
the details of the Gaussian fit to the gas distribution. A 
more common but even less clear analysis uses a dynam- 
ical mass indicator Mdyn ~ AV"^ x R/G, where AV is 
the FWHM linewidth, i? is a radius, and G is the gravita- 
tional constant (e.g. Braine et al. 2001). This quantity has 
the dimensions of mass, but it is not obvious what mass 
it defines, especially when the radius used is not physi- 
cally motivated. However, we can use the quantity to see 
how taking larger or smaller region in the vicinity of ei- 
ther dwarf candidate might affect our conclusions. Specif- 
ically, we evaluate how the ratio of the luminous mass to 
this dynamical mass indicator ch ang es with radius using 
the circular apertures defined in §3.4 and the light curves 
given in Fig. 0. In this case, AV is 2.35 times the one- 
dimensional velocity dispersion, and R is the deconvolved 
aperture radius. The luminous mass is calculated from the 
H I mass and optical luminosity as above. 

The results of this exercise is shown in Fig. |l^. Three 
curves are drawn, each one representing a different stellar 
mass-to-light ratio: the solid curve shows the result for 
M^/Lb = MqLq^, the dotted curve shows the results 
for M^/Lb = 2 MqLq^, and the dashed line shows the 
results for M^/Lb — 5 MqL^. We see that this ratio 
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decreases with larger apertures. Therefore, taking a larger 
radius will not change the conclusions reached above. Tak- 
ing a smaller radius might turn up objects with a larger 
luminous to dynamical mass, but in this case they would 
have a much smaller mass scale than is usually assumed. 

In conclusion, we find no kinematic signature of a dis- 
tinct dynamical entity within the southern tail, but any 
signature might be masked by the strong geometric veloc- 
ity gradients in this region. If there are bound mass con- 
centrations, they require a significant contribution from 
either evolved stars or dark matter to bind them. 

If the lower distance of 13.8 Mpc suggested by Saviane 
et al. (2001) is correct, then the luminous mass decreases 
by a factor of 2 while the virial mass decreases by a factor 
of 1.4, so the luminous to virial mass estimates all decrease 
by a factor of 1.4, making it less likely that either TDG 
candidate is a bound object. Additionally, if there are 
bound regions, the mass scale would be more typical of 
dwarf Spheroidal galaxies than dwarf Irregulars, and the 
gaseous concentrations may have more relevance to rem- 
nant streams around galaxies than the formation of dwarf 
irregulars (Kroupa 1998). 

If there are dwarf-galaxy sized concentrations forming 
within the tail, then the present observations appear in- 
consistent with the tidal dwarf formation scenario sug- 
gested by Elmegreen, Kaufmann & Thomasson (1992). 
In their scenario the velocity dispersion of tidal gas is 
increased over the pre-encounter value due to the inter- 
action, increasing the Jeans mass in the tidal gas and 
leading to the formation of dwarf Irr sized objects. As 
noted in §3.3, the velocity dispersion within the tidal tails 



of NGC 4038/9 is typical of values measured in undis- 
turbed disk galaxies. This is true of long-tailed mergers 
in general (Hibbard & van Gorkom 1996, Hibbard & Yun 
1999a) and agrees with the expectation from numerical 
models. The tails form from the dynamically coldest ma- 
terial on the side of the disk opposite orbital periapses, 
and expand kinematically thereafter. Tails, particularly 
those of low-inclination prograde encounters, should not 
and appear not to be significantly kinematically heated by 
the encounter. And while it is true that an increased ve- 
locity dispersion will suppress the formation of low mass 
objects, it does not follow that it should lead to the col- 
lapse of higher mass objects; it simply means that a higher 
mass is required in order for the concentration to collapse. 
Therefore, the model of tidal dwarf formation suggested 
by Elmegreen, Kaufmann & Thomasson (1992) does not 
appear to apply here. 

4.4. Expected Signature of Merger Remnant 

A last point that we wish to make is that NGC 4038/9 is 
clearly an on-going merger experiencing a massive young 
starburst which is spread throughout the disk of NGC 4038 
and the overlap region. And yet, the long tails of NGC 
4038/9 suggest that the interaction began in earnest of 
order 400 Myr previously. This is very similar to the situ- 
ation in the infrared luminous merger Arp 299 (Hibbard & 
Yun 1999a). Additionally, the proximity of NGC 4038/9 
have allowed the identification of a population of SSCs 
with ages of ^ 500 Myr (Whitmore et al. 1999) spread 
throughout the disk region. Clearly, then, merger induced 
extra-nuclear star-formation began around the time that 
the tails were launched, has continued until today, and will 



likely continue until the merger is complete. Therefore, the 
merger induced "burst" population will be broadly spread 
throughout the remnant, both spatially and temporarily. 
While a significant amount of the present gas supply may 
find its way into the innermost regions of the remnant, 
it is not clear how much of this gas will be converted into 
stars. Much of this gas may be expelled in starburst driven 
superwind (Heckman et al. 1999, Hibbard & Yun 1999b). 

Either way, the expected signature of the evolved merger 
remnant will not be as pronounced as is often assumed 
(e.g., Mihos & Hernquist 1994, Bekki & Shioya 1997, Silva 
& Bothun 1998a, b). Instead of a single epoch metal- 
enriched population confined to the innermost radii of the 
remnant and surrounded by a uniformly-old metal poor 
population, there will be stars with an age range of ~ 1 
Gyr spread throughout remnant body, with perhaps an 
inner population of stars formed from whatever gas is not 
blown out. 



5. CONCLUSIONS 

• We have mapped the Hi in the classical double- 
tailed merger, NGC 4038/9, "The Antennae", with 
unprecedented spatial and velocity resolution. In 
agreement with other observers (van der Hulst 
1979a, Mahoney et al. 1987, Gordon et al. 2001), 
we find the northern tail to be gas-rich at its tip 
but gas-poor at its base, the southern tail to be gas 
rich along its entire length and with a significant en- 
hancement in the gas concentration in the vicinity 
of the TDG. The tail kinematics are broadly con- 
sistent with the interaction geometry inferred from 
numerical simulations. 

• We suggest that the lack of H i at the base of the 
northern tail is due to photoionization by UV pho- 
tons escaping the disk starburst region. If so, the 
ionized gas should be readily detectable, with and 
expected emission measure of order 0.8 cm~^ pc. 

• The H I velocity field at the end of the southern tail 
is dominated by strong velocity gradients which sug- 
gest that at this location the tail is bending away 
from us. The tail velocity gradients may mask the 
kinematic signature of any self-gravitating conden- 
sation in this region. 

• It is not clear whether the TDG candidates identified 
by Schweizer (1978) and Mirabel et al. (1992) are 
self-gravitating. The observed kinematics suggest 
that there is not enough mass in gas alone to account 
for the Hi linewidths of the regions we have delin- 
eated. Further insight requires a measurement of the 
stellar velocity dispersion and tighter constraints on 
the stellar mass-to-light ratio. 

• Whether or not the TDG candidate identified by 
Mirabel et al. (1992) is a self-gravitating entity, it is 
clear that this concentration of gas, stars, and star 
forming regions in the tail is unique. At this loca- 
tion the H I column densities per channel are higher 
than anywhere else in the system, including within 
the main disks. We believe a clue to the origin of 
this concentration is that it occurs just after the the 
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"V" -shape where the two filaments connect back to 
each other occurs just before the location of of the 
TDG. 

• These observations reveal that the "superthin" 
galaxy ESO 572-G045 is a companion to NGC 
4038/9. Our low resolution data show a low column 
density H i extension to the east, suggesting that this 
system may be experiencing ram pressure as it orbits 
NGC 4038/9. 
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Table 1 
VLA Observing Parameters. 



Date 




1996 Jun 3 


1997 Jun 9 


Array 




D 


CnB 


Time on Source (hrs) 




2.5 


6.5 


Phase Center (a, 6, J2000) 




12'' 01" 52'* 


-18° 55' 40" 


Central Velocity (Heliocentric) 




1630 km s-i 




Primary Beam (FWHM) 




30' 




Phase Calibrator 




1130-148 




Flux Calibrator 




3C286 




Correlator Mode 




lA, on-line Hanning 




Bandwidth (MHz) 




3.125 




Number of Channels 




128 




Channel Separation (km s"-'^) 




5.21 




Data set 


High Resolution 


Medium Resolution 


Low Resolution 


"Robust"" parameter 


-1 


+1 


+1, Convolved to 40" 


Synthesized Beam 








— Major Axis x Minor Axis (FWHM) 


11.4" X 7.4" 


20.7" X 15.4" 


40" X 40" 


— Position angle (east of north) 


76° 


24° 


0° 


Noise Level (Icr) 








— Flux Density (mJy beam^^) 


1.3 


0.9 


1.3 


— Column Density (xlO^^cm^^ bcam"-'^ ch~^) 


8.8 


1.6 


0.47 


— Brightness Temperature (K beam^^) 


9.3 


1.7 


0.49 



'^Robust weighting parameter from Briggs 1995. A Robust parameter of -2.5 is equivalent to uniform weighting, while 
a Robust parameter of -1-2.5 is equivalent to Natural weighting. 
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Table 2 

CTIO/UH Observing Parameters. 



Date 


1991 Apr 17/18 


1995 Jan 22 


Telescope 


CTIO 0.9m 


UH 88" 


Camera 


Tek 512 


QUIRC 


Readout Mode 


unbinned 


unbinned 


Focal Ratio 


f/13.5 


f/7.5 


Field of View: 






— Single CCD frame 


3.'8 X 3^8 


3'.2 X 3'.2 


— Final Image 


17'xl7' 


4'x4' 


Pixel Size 


0'.'439 


0'.'189 


Filters 


B, V, R 


K' (2.15 /xm) 


Seeing 


I'.'S 


0'.'6 


Sky Brightness (mag arcsec^^) 


22.1, 21.3, 20.6 




1(7 Sky Noise, binning 1x1 (mag arcsec^^) 


25.5, 25.3, 25.0 




3fT Surface Brightness limit, binning 9x9 (mag arcsec"^) 


27.0, 26.8, 26.5 




Effective Exposure Time 


~2x600s 


3xl20s 



Table 3 

Global Properties of NGC 4038/9. 



Region 


/ Snidv 
(Jy km s~^) 


(Me) 


Velocity Range^ 
(km s-i) 


A R'' 
(kpc) 


Lb 

(Le) 


Mhi/Lb 
(Me L^') 


Disk"* 


>17.2±0.16 


> 1.5 X 10^ 


> 1420-1850 


9 


2.4 X 10^° 


>0.06 


S tail 


32.5±0.23 


2.8 X lO'' 


1590-1770 


65 


3.6 X 10^ 


0.8 


N tail 


4.8±0.10 


4.2 X 10*^ 


1540-1600 


40 


1.3 X 10^ 


0.3 


Total NGC 4038/9 


>54.5±0.49 


> 4.7 X 10^ 


>1420-1850 


110 


2.9 X 10^° 


>0.2 


ESQ 572-G045 


4.3±0.10 


3.7 X 10^ 


1640-1755 


88 


1.9 X 10^ 


1.9 



^Integrated Hi mass, calculated using; M^j = 2.356 x 10"'' J Shi dv Mq, where A is the distance in Mpc, and 
/ Snidv is the integrated Hi emissivity, in Jy km s^^. We adopt 19.2 Mpc as the distance to NGC 4038/9. 

''Range of Hi velocities (heliocentric), taken from the first-moment image. The uncertainty is ±3 km s~^. 
'^Maximum projected distance of the H i from the center of NGC 4038/9. 

'^Observations by Gordon et al. (2001) and Gao et al. (2001) reveal that the cold gas within the disk and overlap region 
emits over the velocity range 1340-1945 km s~^. This is wider than the bandwidth of the present observations, so that 
the disk and total H i fluxes derived here are lower limits. 
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Table 4 

Dynamical Analysis of TDG Candidates 



Parameter 


TDG [MDL92] 


TDG [878] 


Right Ascension" 


12'' 01" 25.7" 


12'' 01" 22. r 


Dechnation" 


-19° 00' 42" 


-19° 00' 12" 


Background" 






Nhi (atoms cm~^) 


3.8x1020 cm-2 


2.9x10^" cm-2 


(mag arcsec"^ 


26.0 


26.3 


Half-Hght Radius", rh (deconvolved) 


3.2 kpc 


3.5 kpc 


Hi velocity dispersion^, a 


13.4 km s-i 


10.5 km s-i 


Mh 


2.4x1OSM0 


1.7xlO«M0 


L% 


8.7x10^Mq 


5.3xlO^M0 




I.IxIO'^Mq 


7.4xlO*M0 


Mi^^/M„r [MJLb - 0] 


0.3 


0.3 


Mlum/Myir [M,/Lb = 2] 


0.4 


0.5 


Mium/M^ir [M,/Lb 5] 


0.7 


0.7 



^Results from fitting a Gaussian plus constant background to the H i integrated intensity map. 
'^Result from fitting a Gaussian to an H i spectrum taken at the peak of the H i distribution. 
"^From the end-points of Fig. |^. See text. 



